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TECHNICAL NOTE 1+J06 


TEMPERATURE AND THERMAL-STRESS DISTRIBUTIONS IN SOME 
STRUCTURAL ELEMENTS HEATED AT A CON3TANT RATE 
By William A. Brooks, Jr. 

SUMMARY 


Analytical, solutions are given for the temperature and thermal- 
stress distributions in thick skins and structural elements such as 
angle, channel, T-, and H-sections when heated at a constant rate. 
Certain of the analytical solutions are evaluated for selected cross- 
sectional proportions. The results are presented in the form of charts 
involving dimensionless temperature, stress, and time parameters and 
therefore are applicable for different materials and heating rates and 
absolute size of the section. The results have been found useful for 
analyzing and correlating experimental data. 


INTRODUCTION 


Interest in the behavior of structural elements exposed to the 
heating rates and elevated temperatures of high-speed flight has resulted 
in a great deal of research, both theoretical and experiment ail, on the 
effects of aerodynamic heating. The variable nature of the heat flux 
encountered in a typical flight plan complicates the problem greatly. 

By considering constant heating rates, it is possible to simplify both 
theoretical and experimental aspects of the problem and still obtain 
valuable information concerning the primary effects of aerodynamic heating. 

The present paper contains one -dimensional solutions for tempera- 
ture and stress distributions in thick skins and structural elements 
such as angle, channel, T-, and H-sections when heated at a constant 
rate. These solutions have been found to be useful for analyzing and 
correlating experimental data. (See ref. 1.) Although some of the 
solutions have been published previously (for example, refs. 2 to k), 
they are derived herein in order to make the presentation complete. 

Certain of the analytical expressions are evaluated for selected 
section proportions. The results are presented in the form of charts 
involving dimensionless temperature, stress, and time parameters and 
therefore are applicable for different materials and absolute size of 
the sections. 




SYMBOLS 


constants 

constants 

specific heat, Btu/lb-°F 
modulus of elasticity, psi 
thermal conductivity, Btu/ft-sec-°F 
length of leg of angle, ft 
integer 

heating rate, Btu/ft 2 -sec 
thickness, ft 
temperature, °F 

maximum temper attire difference, °F 
specific weight, lb/cu ft 
coordinates 

thermal coefficient of expansion, in./in-°F 
root of characteristic equation 
strain, i n./in. 

dimensionless space coordinates 
dimensionless length, Zj./t2 

dimensionless time parameter, kT/cwt 2 
Poisson's ratio 

thermal stress, psi (positive for tension) 
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0 dimensionless temperature parameter, [T - T Q ) — 

V / qt 

Subscripts: 

j junction of elements 

0 initial 

1 element 1 

2 element 2 

av average 

max maximum 

min minimum 

i integer 


TEMPERATURE DISTRIBUTION 


In the following sections, temperature distributions obtained by 
assuming constant material properties and making exact solutions of the 
one -dimensional heat-conduction equation are presented. First, the 
simple case of a skin, thin in the thermal sense, is discussed in order 
to introduce some of the dimensionless parameters which are used in the 
present paper. Then the solutions for thermally thick skins are given 
and the thin skin is examined as a special case of the thick skin. 
Finally, the one-dimensional temperature distribution in some common 
structural shapes is given. 


Thin Skin 

If it is assumed that a skin is heated on one side, that there is 
no variation of temperature through the thickness or in the plane of 
the skin, and that there is no transfer of heat through the unheated 
face by means of convection or radiation, the differential equation 
governing the temperature of the skin can be obtained by equating the 
rate at which the skin absorbs heat to the rate at which heat is provided, 
or 
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, dT 

cwt — = q 
dT 


(1) 


Integration yields 


T - T 0 = 



q dT 


or 


t 3av T 

T - T 0 = 

cwt 


( 2 ) 


where 


‘lav = 7 



q dT 


In a convenient dimensionless form, equation (2) becomes, for a 
constant heating rate. 


0 = e (3) 

where 

0 = (T - To)4 


e - -SL- 

cwt^ 

The principal limitation of equation ( 3 ) is that the variation in 
temperature through the thickness must be negligible; or, in the thermal 
sense, that the skin be thin. Although this solution is called the "thin- 
skin" solution, the actual thickness does not necessarily have to be 
small. A thin skin can best be defined by considering it as a special 
case of a thick plate as is done in the next section. 
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A frequent use is made of equation (l) or (2) in the calibration 
of the heating rate produced at a given point by a heater arrangement. 

The calibration is accomplished by placing a thin skin before the heater 
and, with the heater in operation, recording temperature rise as a func- 
tion of time. Once the rate of temperature rise has been determined, 
the heating rate may be calculated by using equation (l) if the pertinent 
material properties are known. 

Figure 1 shows the relation between the heating rate and the rate 
of temperature rise calculated by equation (l) for some of the more 
common materials. A nomograph of equation (l), presented in figure 2, 
affords a rapid graphical solution of the equation and is applicable 
for more materials than figure 1. The typical material properties 
employed in deriving figure 1 and used in figure 2 are given in table I. 


Thick Skin 

If a thick plate, shown schematically in figure 3> is subjected to 
a constant heating rate on one face, experiences temperature variation 
through the thickness only, and experiences no heat transfer at the 
unheated face, the temperature distribution (derived in appendix A and 
ref. 2) is given by 


4 = e + ^ 


it 2 n=l 


^ 1 ^n cos rutT] e -n 2 ir 2 8 


where 


(4-a) 


0 = (T - * 0 )4 

qt 

8 = -*L_ 

cwfc 2 


T] 



From equation (4-a), temperature distributions have been calculated 
for a thick skin of any material and thickness subjected to a constant 
heating rate. These temperature distributions are plotted in figure 4- 
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as a function of time for given locations through the thickness and in 
figure 5 as a function of location at given times. At large values of 
the time parameter 8(0 > 0.4), the series term in equation (4a) contrib- 
utes a negligible amount to the temperature parameter, which then becomes 
a linear function of time. The temperature-rise rate becomes constant 
for any location through the thickness, and the spacewise temperature 
gradient is independent of time. The result, which can be seen in fig- 
ure 4, is that the temperature distribution becomes fixed, with the tem- 
perature level of the entire slab increasing uniformly. The temperature - 
rise rate at large values of time can also be calculated by considering 
only the mass and heat capacity of the skin and the heat input. 

The boundary conditions of the present problem, namely, constant 
heat input and no loss of heat, preclude the existence of an equilibrium 
condition although a quasi-steady state does exist after an initial tran- 
sient period. The term quasi-steady state is employed herein to describe 
that state in which the temperature-rise rate is constant at all points 
through the thickness of the slab and is hereinafter referred to as steady 
state. In an actual case there are heat losses through both the heated 
and the unheated faces, provided the unheated face does not lie in a 
plane of symmetry or is not insulated, and an equilibrium state does 
exist. The present solution must therefore be considered as an approxi- 
mation for actual cases in which the temperature levels are such that 
excessive heat losses are not involved. 

One disadvantage of equation (4a) Is that it converges slowly at 
small values of 8 . An alternate solution which can be obtained In the 
form of tabulated functions is given by the following equation: 



At small values of 8, the first term of the series suffices. However, 
as 8 becomes larger, more terms are necessary. Although there may be 
some advantage to using equation (4b) for computing temperatures, par- 
ticularly at small values of 8, equation (4a) is more amenable to 
integration and is used to calculate the stresses. 
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The average temperature is found by integrating the dimensionless 
temperature parameter of equation (4a) over the dimensionless thickness 
and the result is 


0av = 9 


■which is the thin-skin solution (eq. 3)* This result is rather obvious 
inasmuch as there are no heat losses and, therefore, all available heat 
is employed in raising the temperature of the mass. If the heating rate 
is constant, then the rate of change of the average temperature must also 
be constant. In figure I|-, the average temperature is shown as a function 
of time by the dashed line. In figure 5 , the spacewise location of the 
point at which the temperature is equal to the average temperature is 
shown by the dashed line. The point at which the temperature is equal 

to the average temperature always occurs in the range ^ tj < 1, this 

ft 

point being at t] = — for the steady state, at which time the tempera- 

' 

ture is a parabolic function of tj. 

The maximum difference in the values of 0 for the outer and inner 
surfaces can be obtained from equation (4a) and is 

4* = | (5) 

From equation ( 5 ) the maximum temperature difference is found to be 


dS! = - S£. (6) 

2 k 

The max imum temperature difference is plotted in figure 6 as a function 
of qt for several materials with properties given in table I. For a 
given heating rate and skin thickness, the temperature difference varies 
inversely as the thermal conduct ivityj that is, materials with the larger 
v alu es of thermal conductivity are associated with the smaller tempera- 
ture gradients. A nomograph of equation (6), presented in figure J, 
permits a rapid graphical determination of the maximum temperature dif- 
ference between the heated outer surface and the insulated inner surface 
of the plate. 
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Another interesting feature of equation (4a) is the time required 
to obtain the maximum temperature difference. Actually, as an examina- 
tion of the series term of equation (4a) will show, the maximum tempera- 
ture difference is reached at an infinite time. Practically, however, 
it can be assumed that the steady state is reached when the series term 
contributes only a small amount to the solution; for example, 2 percent 
of the maximum temperature difference. In figure 8, the time required 
to produce 98 percent of the maximum temperature difference is plotted, 
against the skin thickness. The curves for the different materials are 
parallel straight lines whose intercepts on the time axis vary inversely 
as the diffusivity k/cw. Thus, for a given thickness, the materials 
with the greatest diffusivity require the shortest time to reach the 
condition at which the maximum temperature difference exists . 

One method of defining a thin skin is to set an arbitrary limit on 
the amount of temperature variation that can be tolerated. With this 
method of defining a thin skin, curves such as those shown in figure 9 
can be prepared. In this figure steady-state temperature differences 
through the thickness are given as percentages of the heated- surface 
temperatures. The solid lines represent the difference between the 
temperatures of the heated face and the insulated face. Any combina- 
tion of heated-face temperature T^-^ and qt/k that lies on or above 

the line for the tolerable amount of variation represents a thin-skin 
solution. For example, assume that the temperature of the heated face 
is 1,000° F and that the tolerable amount of temperature difference 
between the two faces is arbitrarily chosen' as 50° F or 3 percent. So 
long as the value of qt/k for the skin under consideration is equal 
to or less than 100, a thin-skin solution will suffice. 


Angle, Channel, T-, and H-Sections 

Consider next an element that is a simplified version of integral 
construction (fig. 10(a)) which may be regarded as part of a skin- 
stringer or skin-web combination. The Idealization employed in the 
present analysis is shown in figure 10(b). All lateral surfaces except 
that being heated are considered to be insulated. Several such elements 
subjected to heating as indicated may be combined to form the sections 
shown in figure 11. 

It is assumed that there is no temperature variation through the 
thicknesses tq and t2 and that there is no heat loss through the 
unheated faces. Because of the discontinuity in the thickness and 
heating, it is convenient to consider the section as two elements as 
shown in figure 10(b), the element exposed to the heat flux q being 
considered as element 2. The following expressions (derived in 
appendix A) give the temperature distributions in the two elements: 
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r 

t 2 

CD 

*1 

El *2 


E 2 tl 




r Si + V 

M *1; 


-Pn 2 ® 


2 _ 

S2 n= 


“ =1 (| + j ) cotp ^-( 1 + lS tmWl 


cos p n (; 
cos 


(7a) 


t l 


0 


Sl + ^2 

ti 


1 ii h .2 i ii ft 2 + 2t 2 
6?^ 2 +25 i 
+ 

51 t 2 

5 2 + *1 


H 


j 2 

t 2 


'Si + u 2 

k S2 


tr 

c 


: 1 2 


,-Pn £ 0 


COS P n Tl 


t 2 S2 n=l /g x t 2 


El *; 


COS 


77 + TtH PnS x - l 1 + ^ d tan 


PnSg 


(7^) 


,52 *1 


where 



(Ti - T 0 ) 


k 

q.t2 



5i- 



*1 


X 

t2 


S 2 “ 


h. 

t2 


4 
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In equations (7) the terms f3 n are roots of the following tran- 
scendental equation: 


to 

sin cos p n ?2 + — sin p n £ 2 cos p n ^ x = 0 

t l 


(8) 


An expression for the average temperature parameter may be found by 
integrating equations (7a) and (7h) as follows: 



(9) 


Upon performing the indicated integration, the following equation is 
obtained: 




il + ^2 

£2 t l 


e 


(10) 


Again, as in the case of the thick skin, the average temperature is a 
function of the heat capacity of the section. 

The infinite series of equations (7) -converge very slowly for small 
values of 9. In order to avoid the tedious procedure of calculating 
the temperature for short times from equations (7)> It is possible to 
obtain closed-form short-time solutions in terms of tabulated functions. 
The following equations (derived in appendix A) are approximate short- 
time solutions for equations (7a) and (78), respectively: 



(11a) 
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02 ■ 0 - . ~ <j | 6 + 


(£2 “ Tl)\||- ex P 


(£2 - ^ 


erfc 


^2 ~ T 

2^8 


(£2 - *0 

14-0 


-]4 e + ^ 


tl)' 


erfc! 


£2 + *1 > 

2\[e 


(£ 2 + ti\Jz exp 


(£2 + 7 l) < 

46 


(Hb) 


Although these short-time equations give good approximations to the 
maximum fmd mini mum temperatures , they do not satisfy continuity at the 
junction of the two elements. For short-time solutions a satisfactory 
approximation for the temperature at the junction is 





( 12 ) 


Of course, for long times the series terms in equations (7) become 
negligible and a quasi-steady state is reached in which the temperature 
is a lin ear function of the time parameter. The slope of the linear 
relation between the time and temperature parameters is, at steady state, 
determined by the heat capacity of the section. The existence of a 
state of equilibrium is precluded by the initial, assumptions as was the 
case for the thick plate. 

The maximum temperature (at tj = 0), the minimum temperature 
(at £ = 0), and the temperature at the junction of the two elements 
(at £ = £ 1 , T] = £ 2 ) have been calculated for twenty-seven cases (all 

possible combinations of — — = 1, 2, and 4) £ 2 . = 5? 10? and 20} and 

t x 
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^2 = 10, 20, and 30 j. In figure 12 the maximum temperature is plotted as 
a dimensionless temperature parameter 02 max a gs-i&st the dimensionless 

time parameter 6, each part of the figure, being for a different thick- 
ness ratio t2/tp. In figure 12 a straight line with slope of 1 is identi- 
fied as = 0. This line is the same as the "thin-skin" solution and 

serves as the upper limit for the maximum temperature by representing the 
case in which the maximum temperature of the heated element is not influ- 
enced by the presence of the unheated element. This line can also be 
represented as £2 = °o, but in either case, £]_ = 0 or £2 = “, the 

thickness t£ cannot he zero or large enough to violate the assumption 
that there is no gradient through the thickness t2- If £2 suffi- 
ciently large, say of the order of 50, the maximum temperature of the angle 
section is only slightly affected by the presence of the unheated leg. 

In figure 13 the temperature at the junction of the two elements is 
plotted as a dimensionless parameter 0j against the time parameter 0. 

If Z-^ = Z2, the temperature parameter f j is a linear function of 0 

given by the following relation: 


Hr 6 


(15) 


1 + 


which is a special case of the average temperature parameter 


tg 

0 = — 0 

rav ^ 

ja + ^2 

£2 t i 


(14) 


Zp 

For those cases in which — — > 1, the lower limit is the straight 


line defined by equation ( 13) - For those cases where — < 1, the upper 

*1 

limit is the line defined by equation (l3)» 
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In figure l4 the minimum temperature is plotted as a dimensionless 
temperature parameter against the time parameter 6. The upper limit is 
the thin-skin solution which corresponds to = 0., or = 0, and the 

lower limit is the horizontal axis which corresponds to = «. The 

minimum temperature of the angle section is more sensitive to changes 
in the value of than to changes in the value of £ 2 * 


THERMAL STRESSES 


In the following discussion, it Is assumed that the section being 
considered is sufficiently far from the ends of the structural element 
that conditions at the ends do not affect the stress distribution. The 
stresses are determined by elementary elastic theory employing the 
assumptions that plane sections remain plane and that material prop- 
erties do not change with temperature. The derivations of the stress 
equations are given in appendix B. 


% 


Thick Skin 

The stress distribution in an unrestrained thick skin is given by 
the equation (see appendix B and ref. 4) 


a = 


otE 


1 - H 


(T 


t 

- T°) - i J* (T - T 0 )dy - 12 i ^ (T - T Q )y dy 


( 15 a) 


where 


T = T(y) 


with y being measured from the midplane of the plate. If y is meas- 
ured from the unheated face of the plate, equation ( 15 a) must be replaced 

*>y 


a = - 


aE 

1 - H 


(T - T 0 ) + (T - T 0 )dy + - A A (T - T 0 )y dy 

t\t 3/^0 t 2 \ 2 v J 0 


( 15 b) 
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Symmetrical heating .- If the skin is heated at the same rate on both 
surfaces, equation (ipa) reduces to 


0 


ctE 

1 - H 


(T " T 0 ) 


2 

t 



( 16 ) 


When equation (4a) is substituted into equation (l6) and the indicated 
integration is performed, the result is 


l— = _ j£ + i + if" cos gtij e -n2jt2e 

/qt 2 6 jj 2 n 2 


(17) 


Stress distributions calculated by equation ( 17 ) are plotted in fig- 
ures 15 and 16 in terms of dimensionless parameters; and, therefore, 
the results are applicable to any material for any thickness or time. 

In figures 15 and 16, T] = 0 corresponds to the midplane of the plate. 

The stress distribution is symmetrical with respect to the midplane; 
therefore, the stresses are shown only for the half -depth of the plate. 

Asymmetrical heating .- If the skin is heated on one surface only, 
the resulting unsyrametrical state requires that all terms of equation ( 15 b) 
be included. When the temperature-distribution equation (4a) for this 
condition is substituted into equation (l^b) , the resulting stress equa- 
tion is 


0 



tl(l - tj) - 



(_2jn cos DJtt) e -n 2 rt 2 0 


n c 


+ 


ifd - 2n) It 

n=l,3,5, 


.-n2jt20 


nr 


( 18 ) 


Stress distributions calculated by equation (l8) are plotted in 
figures 17 and 18 in terms of dimensionless parameters. It Is interesting 
to note that, in certain regions of the skin, the thermal stresses reach 
maximum values at relatively short times during transient heating. In 
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other regions the maximum stress is not reached until steady state exists. 
The inclusion of the last term in equation (l^b), which term, is obtained 
by satisfying moment equilibrium, causes the stresses to reach maximum 
values at relatively short times. 


Angle, Channel, T-, and H-Sections 

If it is assumed that material properties do not change with tem- 
perature and if elementary theory is employed, simple expressions may be 
derived for the stresses acting on unrestrained sections such as those 
shown in figure 11 with the heating as indicated. The stress equations 
are in terms of geometric properties, integrals involving the tempera- 
ture distribution, and the local temperatures and are as follows: For 

the H-section, 


£1 = 
ctE qt 2 




( 19 a) 


ccE qt 2 



0 2 


(19b) 


for the T-section, 


fl k 
otE qt 2 




+ iJLk + *M± 
1^2 *ij\ h 



h 


(20a) 



1 6 
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^2 k 
aE qt 2 


6 2 




*1 


6i 



" 02 


for the channel section, 


fl k 
aE qt 2 



t l $2/ 



z' 


f2 k 

otE qt 2 





II + 



(20b) 


(21a) 


(21b) 


for the angle section, 


fl Jl_ 

aE qt2 





2 Cl(— - -Vl + 

2(2 - 1\ + — M/jl _ ij 

\h v 

_ \ £1 / t i ^i\fi / 


V 


—(3 -L 

S 2 V «i 



- 0i 




(22a) 
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IT 


S2 k 
ccE qt 2 


k . *2 


+ 


£0 t. 


2 ^(& ' 3) 11 + ' 5 ^ " ^ 12 + 


52 

t i k 


22 


JL - l\ + !i 2i(x . 

k J t 2 ^2\j2 


z 3r “02 


The integrals I;j_, I 2 , and I 5 are defined as follows: 


(22b) 


*i - 


qt 2 


(23) 


where 


Cl' = r 1 (T l - T Q )d£ + ^ r' 2 (T 2 - T Q )dTi 

Jo t x Jo 


V = P < T i - T o)U - 5i)d5 

dO 

£ 

I 3 ' = f 2 (t 2 - ToJCn - 5 2 )4 ti 

Jo 


-/ 


(24) 


The quantities I 2 and Ij have been calculated for all possible 


combinations of 


— = 1 , 2, and 4, = 5 * 10 > and. 20 , and £ 2 - 10 , 20 , 

t l 

and 30 ) and are plotted in figures 19 and 20. Because of its simplicity, 
Ij_ is given only in the form of an equation. (See appendix B.) 
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EXAMPLE 


Temperature Distribution in an Angle 


In order 
an example is 


to illustrate the results of the analysis described herein, 

presented. The case chosen is ^2. =1, £1 = £2 = 10* For 

*1 


the example chosen 
eq. (8)) 



and the characteristic values are (from 


3n = — ( n = 1> 3j 5; • • • ) 

2 ^1 


The temperature distribution was calculated by using equations (7) and 
is given in figures 21 and 22 as a function of time and space, respec- 
tively. This particular case reaches quasi- steady- state conditions 
rather rapidly (0 » 150). Note also that the temperature at the junc- 
tion of the two elements that make up the basic angle section Is a linear 
function of time and is the average temperature of the section. 


Stress Distribution 

The quasi -steady-state stresses were calculated for the four cases 
shown In figure 11 by assuming 0 = 00 and using equations (19) to (22). 
The results are presented in dimensionless form In figure 23* The 
results, which are for unrestrained expansion. Indicate that, for the 
case chosen, if the H-section is divided along either line of symmetry, 
the absolute value of the maximum stress on the resulting section Is 
one-half that of the H-section. If the H-section Is quartered, the 
absolute value of the maximum stress on the resulting section is one- 
fourth that of the H-section. 


CONCLUDING REMARKS 


Solutions for the temperature distributions and thermal stresses of 
structural elements such as plates, channel, angle, T-, and H-sections 
when subjected to constant heating rates have been presented. These 
solutions and charts obtained from them have been found to be useful in 
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analyzing and correlating experimental data. In addition to temperature 
distributions, thermal stresses are given in forms involving simple 
integral functions of the temperatures. 

The temperatures of selected points on the element cross sections 
have been calculated and presented in dimensionless form for several 
geometric configurations. The temperature integrals involved in the 
thermal stress equations have been calculated for the same configurations 
and are also presented in dimensionless forms. By the use of dimension- 
less temperature and stress parameters it is possible, to employ the 
results for different materials and heating rates . 

A detailed distribution of temperature and stress is given for a 
selected case in order to provide some insight into the nature of the 
variation of temperature and stress on the cross section. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 27, 1958 . 
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APPENDIX A 


TEMPERATURE DISTRIBUTION 


Thick Skin 

In figure 3> a schematic diagram of a thick skin heated on one face 
is shown. It is assumed that there is a variation of temperature through 
the thickness but no variation in the plane of the skin. The governing 
differential equation is 


ST _ k S 2 t 
St cv Sy 2 


(Al) 


In addition, it is assumed that there is a constant heat flux on 
one face (y = t) and that there is no heat loss from the other face 
(y = 0). The resulting boundary conditions are 


ST 

dy 


t,T 



(A2a) 


$1 . 0 

^y o, T 


(A2b) 


The initial condition is 


T(y,o) = T 0 


Expressed in a dimensionless form, equation (Al) becomes 


(A3) 


St = S^t 

Se " s n 2 


(A4) 
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and equations (A2) are "written as 


St 


1,0 



(A5a) 


M = 0 
^ 0,0 


(A5b) 


where 


0 = J£_ 

cwt2 



(A6) 


By employing the Laplace transformation the differential equa- 
tion (A4) and the boundary-condition equations (A5) may be written as 


and 


— - sT + T 0 = 0 
dT] 2 


(A7) 


dT = qt 
dT] ■]_ 


(ABa) 


dT 

dq 


= 0 


(ABb) 
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where 


t(tj,s) = l {r(t|,e)} 


The solution of equation (A 7 ) subject to the boundary conditions 
(A 8 ) is 


T = ^ 2 . + S& cosh \/bt) 
5 k s 5 / 2 sinh 


An inverse transform of equation (A10) and the required temperature 
expression is 


T = T o + ^S.0 + H--i - — (_i)n cos ratt) ^n 2 ^ 2 © 


2 6 « 2 n=l 


An alternate form of solution which converges rapidly for small 
values of 0 is obtained by writing equation (AlO) in the following 
equivalent form 


t = 5a + ai 1 

S k s 5 / 2 n=l,3,5, . . 


[ e -\/s(n-Ti) + e -^/F(n+Tj)J 


The inverse transform of equation (A12) is 


T = T~ + — 




n=l,3,5 } 


(n - T])erfc[ n ~ T l | + (n + T])erfc(- 

\2\JQ I \2 \fd 


r 

fe + 6X5 - Ie l)2 . 


9 
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The average temperature parameter of the thick skin is found by 
rearranging equation (All) and integrating over the thickness, or 


yielding 



0 av = 6 


(AllO 


(A15) 


The nondimensional temperature difference, defined as 

4 - f*|i,e ' 0|o,e 

may he found by using equation (All) and is 


± (ALT) 

2 *2 n=l,3,5, ... n 2 


It is obvious that, as 0 becomes large, the nondimensional tem- 
perature difference increases until it reaches its maximum value of 


A0 = \ (A18) 

or 

iff = (A 19 ) 

2 k 


at an infinite time. 
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Angle, Channel, T-, and H-Sections 


The idealized angle section of figure 10(b) is next considered. 

It is assumed that there is no variation of temperature through the 
thicknesses t^ and t 2 and that there is no heat loss through the 

unheated faces. Because of the discontinuity in heating and thickness, 
the section is considered as two elements wherein the temperatures are 
governed by the following equations: 

_ (k) ^1 

W/j. ay 2 

(A20a) 

&S U \ a%2 , 4 

3t Wa *2 («rt) 2 

(A20b) 

subject to the following boundary conditions 


dr,. 

— “(0,T ) - 0 
dy 

(A21a) 

dr, St p 

Ort) x = -(tb) 2 — ~{l2>r) 

ay ax 

(A21b) 

T l(il^r) = T 2 ( Z 2 # t ) 

(A21c) 

dPo . . 

— g (0,r) = 0 

Sx 

(A21d) 


The initial condition is 


T L (y,0) = T 2 (x,0) = T q 


(A22) 
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If it is assumed that the two elements are of the same material, 
the following dimensionless equations may he written to replace equa- 
tions (A20) and (A2l): 


St 1 

S 0 ££2 

&T 2 5^12 Q.^2 

Se 2 k 


and 


^(0,6) . 0 

35 


| - - 


ar 2 

an 


(C2,e) 


T L (5i,e) = t 2 (£2, 0 ) 


where 


% 



J 


(A23a) 

(A25h) 

(A24a) 

(A2lfb) 

(A24c) 

(A24d) 

(A25) 
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If the Laplace transformation is employed, the differential equa- 
tions (A23) and the boundary-condition equations (A2b) become 


+ 

t-3 

O 

II 

O 

(A26a) 

d^Tp — 1 qto 

- sTo + T 0 + ± = 0 

j 2 ^ 0 s k 

d n 

(A26b) 

0 

11 

0 

(A27a) 

*1 dT l,, x dT 2 f, , 

— — (£ 1 ) = ■ ~^2) 

1 2 d£ dr} 

(A27b) 

\(ii) = t 2 (c 2 ) 

(A27c) 

dT 2 N 

— — ( 0 ) = 0 

d *l 

(A27d) 


The solutions of the transformed differential equations (A26) sub- 
ject to the boundary conditions (A27) are 



qt 2 

s^k 



to 

sinh yfs^ cosh + — sinh cosl1 


(A28a) 



NACA TN 3^06 


27 



sirih Vs£i cosh \fir) 


sinh ^s5 1 cosh y[s£ 2 + “ si 11 * 1 2 cos ^ L \f®^l 


(A28b) 


The poles of equations (A28) are a pole of order 2 at s = 0 and 
simple poles at s = -p^, where the terms P n are roots of the following 
transcendental equation: 


■fcp 

sin P^p cos P^ + £j- sin Pn^2 cos Pn^l = 0 


(A29) 


For cos cos Pn(; 2 ^ Pn is given by the following equation: 


tan p^q + ^ tan P n £ 2 = 0 (A30) 

il 

In the present paper the values = 1 and 3 occur. For the case 

5a. 

£ 2 

— = 1, P n is given by 

^1 

cos PnCi =0 (n = 1, 3 > 5* • • •) 

or 



(n = 1, 3, 5, • • •) (A31) 
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£ 2 

For — = 3> there are three sets of p n . Two of the sets are 

^1 

roots of equations (A30) and (A3l), and the third set is given by 

cos = 0 (n = 1^ 5* 7 > 11, • • •) 

or 


Pn - 


rut 

6 Sl 


(n = l, 5 , 7, 11, ... ) 

(A32) 


By a formal application of the inversion integral the required tem- 
perature distributions are found to be 


(T X - T 0 ) 


qt 2 


^2 

*1 


e 


1 ¥ 

2 * 




Cl + H Cl + t2 

C 2 t i C 2 t-L 


'Ci + 

■ C 2 Hi 


^ 22 - 


,-Pn 2 6 

Pn 5 


cos p n f; 


f-2- 7 

^2 n=l /£, t 2 \ / t 2 \ cos p n ^ x 

'r + r- cot - 1 + t 1 r taD p n?i 


C 2 Hi 


Cp 


(A33a) 
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(T 2 - T 0 )-S- - ^ 


e 




£1 J. ±2 

£2 *1 


5l + ^2 

^2 *1 


ti + i 2 

^2 t ] 


t l 2 




fZ- 


.-Pn 20 

3n 5 




|tan p n ^ 2 


cos P n T] 
cos ^2 


(A35b) 


Because equations (A33) converge slowly for short times, it is 
advisable to derive a short-time solution for the present problem. This 
derivation can be made most readily by appropriate substitutions in equa- 
tions (A28) . For small times (large values of s) assume that 


and 


sinh n/s - ^ = cosh \fs£-j_ = ^ e^^l 

> 

sinh \j"s£ 2 = c °sk \[^2 = ^ e^^2 


cosh \Jb£ 


cosh \fsT] 


l( e #? + e -\P) 

► 




y 


(A34) 


(A35) 
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If the q uan tities defined in equations (A35) are substituted into 
equations (A28), the result is — 


*1 


To 

s 


qt 2 

k 


i + h 

t2 


_1_ 

c 2 





t 2 




(A36a) 


(A36b) 


The inverse transforms of equations ( A 36 ) can be written as (ref. 5) 


(T X - T 0 )-J- = 
qt 2 


1 + 


t 2 


9 + 


(h - if 


erfc 


^ 

V 2 \/~Q 1 



(A37a) 
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(t 2 - t 0 )J£- = e 
qt2 




(A37b) 
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APPENDIX B 
THERMAL STRESSES 


In the following discussion, secondary effects are neglected; that 
is, it is assumed that the section under consideration is sufficiently 
■removed from the ends of the structural element that the stress distri- 
bution is not influenced by the stress-free condition at the ends. All 
the stresses are determined by elementary elastic theory employing the 
assumptions that plane sections remain plane and that material properties 
do not change with temperature. 

If plane sections remain plane, the axial strain at any point can be 
written as 


e - a + bx + cy 
The stress-strain relation is 


(Bl) 


6 = £- + o(T - T 0 ) (B2) 

E* 

where, for a free plate with temperature variation through the thickness 
only, 

E* = E (B3) 

1 - p 

and for plane stress 


E* = E (B4) 


The stress can be expressed as 


8 


a = E* a + bx + cy - a(T - T ( 


>>] 


(B5) 
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The equilibrium equations which must be satisfied are 



(B 6) 


If equation (B5) is substituted into equations (B 6) and the indi- 
cated integrations are performed, there is obtained the following general 
set of equations for determining the values of the constants a, b, 
and c : 


A^a + B]b + C-^c 
A^a- + B 2 b + C 2 c 
A^a + Bjb + C^c 


a 'I 1 ' 


aI 2 ' > 

al 5 r 


(BT) 


The constants A, B, and C are known functions of the geometry 
of the particular section under consideration; I^’ , I 2 ’ , and 1^' are 

integrals involving the temperature . 


Thick Skin 


The thick skin considered in the present paper experiences varia- 
tion of temperature in the thickness direction only. The coefficients 
for equation (B7) for the free plate are 


A! = 1 
a 2 = 0 

a 5 = i 
0 2 


B]_ = 0 
B 2 = 0 

b 5 = 0 


Cl 

c 2 

c 


t 

2 

0 


3 = T“ 


t_2 

3 


(B8) 
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The right-hand sides of the equations (B7) are 


aI i’ = r /■* ( T -' T o)dy 
t Jo 


al 2 ' = 0 


al*’ = & (T - T 0 )y 

t J 0 


dy 


(B 9 ) 




When equations (B7) are solved for the coefficients a and c 
which are then substituted into equation (B5) , the stress is 


cr = 


ocE 


1 - fi 


(T - T °> + IH 


)fo (T - 


T 0 )ay + 


Si - i)fo (t - T ° ,y * 


(BIO) 


where y is measured from the unheated surface . If the origin i6 trans- 
ferred to the midplane of the plate j equation (15a) Is obtained. 

If equation (All) is substituted In equation (BIO), the resulting 
stress equation for the unsymmetrically heated plate Is 


LdJiW (1 ) _ 1 + 4_ (-D n —-W e“ n2lt20 + 

^ ctE y qt 6 ^2^1 n 2 


- 2rj) tl 

^ n=l, 3 i 5 ,. 


n^ “ 


(Bll) 


If the skin is heated symmetrically with respect to its midplane 
B.nfi there Is no variation of temperature in the plane of the skin, the 
last two of equations (B 6) are automatically satisfied and b = c = 0. 
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The first of equations (B6) must "be employed to determine a. When 
this is done, the stress is 


a 


aE 

1 - H 


(T - To) - 7 

"C 



(B12) 


where y is measured from the midplane of the plate. 


If equation (All) is substituted into equation (B12), the resulting 
stress equation for the symmetrically heated skin is 

af 1 - As . . 2 (-I)B <=°a Ml e -n 2«29 , } 

\ aE / qt 2 6 n a 


Angle, Channel, T-, and H-Sections 


The thicknesses of the sections of figure 11 are small enough to 
permit the assumption of plane stress. For the general case, the coef- 
ficients A, B, and C of equation (B 7 ) are as follows: 


ai = Safe ♦ 

Bl 


■s 

Oi - ^ + ^ 

V s V 

\ft 2 \j 

2 ^ ^2 \) 

q 2 

Ae = - — 

b 2 

.. qq 5 

A 

CM 

CM 

H 

CM 

+3 

1 

II 

CM 

O 

2 

6 

2 

t 2 

A ^2_ 

b 3 

t 2^1^2 2 

„ t 2^2 5 

5 2 

2 

c 3 .- 6 

J 


(B14) 
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The right-hand sides of equations (B7) are 




aI l' » a P 1 (T 1 - T 0 )d£ + ^2 r^ 2 (t 2 - T 0 )drj 

y 0 ti jp 


r 

al 2 ' = a / (T X - T 0 )(g - ^)d^ 
J 0 


£ 

alj' = a T 2 (T2 - T 0 )(r} - 

Jo 


J 


(B15) 


The H-section of figure 11 is the simplest case because of double 
symmetry. In this case the constants b and c are zero; therefore. 


£l 

aE 


Mr + - 
V=2 *1, 


- (Tl - To) 


(Bl6a) 


£2 

aE 


II' 


♦ * 
V 2 *i ; 


- (T 2 “ T 0 ) 


(Bl6b) 


For the T-section of figure 11, the constant b equals zero and 
the stresses are 


£l 

aE 


MS -S’ 


_L 

Ui 


V. 


it 


>1 *- 


'*JL + V 


V^2 


12 i- 


I h ' + 


t 2 

11 -r 

t i 


v 


j 


- (Tl - T 0 ) 


(B17a) 
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£2 

oE 


/£i , t2\\3 x e n 

^ 2 L + * 


- (t 2 - T o) 


For the c hann el section of figure 11, the constant c 
and the stresses are 


£l 

aE 


- In ’ + — 


t 2 I 


£]_ ^ 2 \ l ^ ^1 ^2 


- (Tl - T 0 > 


t 

*2 *1, 


£2 

aE 


6 


K 4 I + S 

i - 1\ 
s 2 * 1^2 ! 


ii 


JL _ i> 

fip 3; 


II* + 


2 ^i_/jL _ 1 1 h 

£2^ 


l3’ ? - ( T 2 " ®o> 


The angle section of figure 11 requires the use of all 
and the stresses are 


£l 

aE 


!A (H, 

£2 Wl_ _ li 

* 1 S 1 U 1 > 


- tV + 


2f 2 i- - 1 \ + 
£ 1 




»' + — —(3 — - iW 


tl^l £l 


(B17b) 


equals zero 


(Bl8a) 


(Bl8b) 


the constants. 


- (Ti - T 0 ) (B19a) 
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r 


£2 

aE 


V 2 l^ + - 






+ 




+ h k/a . 

t 2 £2^2 



(T2 - T 0 ) 


(B19b) 


In order to keep the results in a dimensionless form, rather than 
employing the integrals defined by equations (l5)> integrals defined as 
follows are used: 


li = Ii‘ -7- 
qt 2 


(B20) 


After the indicated integrations are performed, the equations for 1^, 
I 2 > and I3 are as follows: 


to 

ii = r 

t i 


(B21a) 


I 2 - - 


*2 hL 

tl 2 


6l + ^2 
^2 t l 


e + 


12 


1 2 + 2 f 2 ) + 1 r 2 ts 

s + 25 2 ; + 5 ti ^ 


Si t 2 

^2 t l 


,-3n 2 9 


Pn 5 


^2 

^1 ^2 H-l / to\ ( tp\ 

^ + -7 «* M2 - 1 + r r tan 
^2 t l/ V ^2 *1/ 


'1 - 


cos Ml; 


(B21b) 
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The stresses may now he written in a dimensionless form by multi- 
plying the left-hand sides of equations (Bl6), (B17), (Bl8), and (B19) 
by k/qt 2 and replacing I^ 1 in the right-hand sides by Ij_. The 

"bp 

values of Ip and It: were calculated for all combinations of — = -l, 

5 ti 

2, and 4, £-j_ = 5 , 10, and 20, and £ 2 = 10* 20, and 3 0 , by using either 

equation (A 30 ), equation (A 31 ), or equations (A30), (A3l), and (A32), 
which define the characteristic value j3 n . For a given case 1^ can 
be computed from equation (B21a) . 
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TABLE I.- TYPICAL MATERIAL PROPERTIES 


Material 

T, °F 
(a) 

w, lb /cu ft 

c } Btu/lb-°F 

k, Btu/sec-ft-°F 

Aluminum alloy 

250 

173 

0.220 

0.0228 

Magnesium alloy 

250 

110 

ro 

.0145 

Copper 

4 oo 

570 

.097 

.0583 

Titanium alloy 

500 

280 

.145 

.0021 

Stainless steel 

600 

490 

.130 

.0031 

Monel 

700 

520 

.118 

.0044 

Inconel X 

800 

505 

.128 

.0028 


Representative of average temperatures in typical applications. 
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Figure X.- Relation between heating rate and rate of temperature rise 
for thin plates of various materials.' (See table I for material 
properties . ) 
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Btu/f t^-°F ft 


Ijioo J 
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-- 
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Inconel x-^^0-- 

-‘-l 

, 1000-1. 


+.001 


dT 

Key: The lines 1 and q-cv intersect on P-P. 

dT 

Figure 2.- Nomograph for equation relating heating rate and rate of 

temperature rise in thin plate, q = cvt — . 

dT 
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Figure 4 .- Temperature history of thick plate subjected to constant 

heating rate. 
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Figure 5*~ Temperature distribution within thick plate subjected to 

constant heating rate. 
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Copper- 


Aluminum alloy- 
Titanium alloy- 


Btu/ft-sec 


-Magnesium alloy 


Monel- 


Stainless steel- 
inconel X 


50 100 

AT, °F 


Figure 6 .- Effect of constant heating rate on maximum temperature dif 
ference between heated and unheated surfaces of thick plates of 
various materials. (See table I for material properties.) 
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2 AT , q, 

°F Btu/ft.^-sec 

1000-rl00 


100 


10 -LI 


10 


k, t, 
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.1 T 1 
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.01 -j-.l 


Monel — 
Stainless steel 
Inconel X' 
Titanium alloy — I 


.001 - L .01 


Key: The lines 2Z£F-k and q.-t intersect on P-P. 

Figure 7 .- Nomograph for equation relating heating rate and maximum 

temperature difference through thick plate. a 3 : 5 i& 

2 k ‘ 
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t , in . 


Figure 8.- Time required to produce 98 percent of maximum temperature 
diff^rence between heated and unheated- surfaces of thick plates of 
various materials subjected to constant heating rate. (See table I 
for material properties.) 
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Figure 9 .- Assigned temperature differences as function of temperature 
of heated surface of plate and heating rate. 
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Figure 11.-. Structural shapes formed by combinations of angle sections. 
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Figure 17.- StresB distribution in thick plate subjected to constant 

heating rate on one surface. 
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Figure 19.- "Values of function Ig appearing in stress equations 
for angle, channel, and T-sections. 
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Figure 19*- Continued. 
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Figure 19 . - Concluded. 
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Figure 20.- Concluded. 
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